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Abstract: We study the constraints imposed by neutrino oscillation experiments on the 
minimal extensions of the Standard Model (SM) with ur gauge singlet fermions ("right- 
handed neutrinos"), that can account for neutrino masses. We consider the most general 
coupling of the new fields to the SM fields, in particular those that break lepton number 
and we do not assume any a priori hierarchy in the mass parameters. We proceed to 
analyze these models starting from the lowest level of complexity, defined by the number of 
extra fermionic degrees of freedom. The simplest choice that has enough free parameters 
in principle (i.e. two mass differences and two angles) to explain the confirmed solar 
and atmospheric oscillations corresponds to tir = 1. This minimal choice is shown to be 
excluded by data. The next-to-minimal choice corresponds to tir = 2. We perform a 
systematic study of the full parameter space in the limit of degenerate Mai or ana masses 
by requiring that at least two neutrino mass differences correspond to those established 
by solar and atmospheric oscillations. We identify several types of spectra that can fit 
long-baseline reactor and accelerator neutrino oscillation data, but fail in explaining solar 
and/or atmospheric data. The only two solutions that survive are the expected seesaw and 
quasi-Dirac regions, for which we set lower and upper bounds respectively on the Majorana 
mass scale. Solar data from neutral current measurements provide essential information to 
constrain the quasi-Dirac region. The possibility to accommodate the LSND/MiniBoone 
and reactor anomalies, and the implications for neutrinoless double-beta decay and tritium 
beta decay are briefly discussed. 
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1. Introduction 

The existence of tiny neutrino masses is probably the first signal of physics beyond the 
Standard Model (SM). One can envision various possibilities for such new physics, but 
the simplest is probably the addition of extra singlet fermions, that can play the role 
of the missing helicity partners of the SM neutrinos. Such models include very different 
possibilities ranging from Dirac neutrinos to the popular Type-I seesaw models j|, pj, ||, |4j . 
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The particle content is the same in both cases, but the couplings of the extra especies to 
the SM fields and among themselves are different. 

In the case of seesaw models, the most general couplings are assumed, in particular 
those that break lepton number, such as Majorana masses for the singlet fermions. In order 
to recover naturally light neutrinos, those masses are assumed to be very large, above the 
electroweak scale. There is however no reason why this scale should be so high, and 
several works have considered the phenomenological implications of a much lower seesaw 
scale H ll 0, g. Dirac neutrinos can be recovered by imposing a global lepton number 
symmetry. Such a lepton number symmetry can however be implemented in more than one 
way, since there is freedom in the lepton number charge assignments of the extra sterile 
species. Different choices of these charges led to different models besides Dirac neutrinos, 
such as the inverse seesaw H [R| . 

The purpose of this paper is a first step towards a systematic exploration of the phe- 
nomenology of such models, in increasing order of complexity, in order to quantify the 
constraints imposed by data. The theoretical prejudice will be just the number of extra 
Weyl fermions, ur. We consider therefore the more general couplings, i.e. without im- 
posing lepton number symmetries, although obviously the models with exact symmetries 
can be recovered in appropriate limits. In this work we will consider just the simplest 
possibilities hr = 1, 2 and leave other cases for future investigations. 

The simplest choice, which corresponds to adding just one additional singlet Weyl 
fermion, ur = 1, has in principle enough free parameters to fit the solar and atmospheric 
oscillation data. A more detailed analysis however shows that it does not work, just 
considering long-baseline oscillation neutrino data. The next-to-minimal choice requires 
two Weyl fermions, ur = 2. Such a possibility is of course well known to be a viable 
one (as long as the LSND/MinBoone signal is discarded), both in the Dirac limit and 
in the standard seesaw limit. In both cases, the physics spectrum contains one massless 
neutrino and two massive ones and, therefore, the two required splittings. The main goal 
of the paper will be the exploration of the parameter space in between these two limits, to 
search for other viable solutions that could accommodate at least the solar and atmospheric 
oscillation. In order to analyse this model we discuss a convenient parametrization valid 
for ur < 3. In the mini-seesaw region, the Casas-Ibarra parametrization [11] is used to 
derive approximate oscilllation probabilities, accurate in this regime. 

Many works before have been devoted to study the implications of neutrino oscillation 
data on models with extra sterile neutrinos (some recent analyses are [12, 13, 14, jig]), 
usually refered to as 3 + 1, 3 + 2, ...3 + N s . Most of these studies have been done with 
the motivation of trying to accommodate LSND |[(J, and MiniBooNe data [|l7], [Djf| . It is 
important to stress that these phenomenological models usually correspond to a generic 
model with 3+N s mass eigenstates. The number of free parameters for N s = ur is typically 
much larger than what we find, either because the number of Weyl fermions involved is 
different (e.g. 3 + 1 Dirac fermions correspond in our context to ur = 5 and not to ur = 1) 
or because couplings that are forbidden by gauge invariance in our model, such as Majorana 
mass entries for the active neutrinos, are included effectively in the phenomenological 
models. It is obvious that our models for any hr will be contained as restricted cases in 
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these phenomenological models for N s = tlr. However, the analyses performed in these 
works do not take such limits and usually restrict the number of parameters by assuming 
instead some hierarchies between neutrino masses that could accommodate LSND. 

The structure of the paper is as follows. We start by setting up the notation in sec. || 
and the parametrization in sec. [3[ In sec. |3| we consider the implication of oscillation 
experiments on the 3+1 model. In sec. |5| we consider the 3+2 model in the degenerate 
limit in the context of LBL, solar and atmospheric oscillations. In sec. || we discuss the 
implications for LSND/MiniBoone, as well as for tritium and neutrinoless double beta 
decay experiments. In sec. we briefly consider the non-degenerate case in the mini- 
seesaw regime. The very stringent bound on the quasi-Dirac region implied by solar data 
is discussed in detail in sec. |8[ 

2. Models with extra singlet fermions 

The leptonic part of the Standard Model Lagrangian is: 



a=l 



li-fO^L + ^D^el - £ liX^4 + h.c, \ (2.1) 

/?=! 



where l a is the lepton doublet. To simplify notation, from now on summation over repeated 
indices is understood. 

We study the minimal extension of the Standard Model that can account for neutrino 
masses, i.e. the inclusion of tir gauge singlet Weyl fermions (hereafter called "right-handed 
neutrinos", z/r): 

C = C SM - ItYf^n + h.c. (2.2) 

If lepton number is violated, gauge-invariant Majorana mass terms can be added to the 
minimal Lagrangian: 

C = C SM - ItYf^n - \vtMfv^ R + h.c. (2.3) 

The new mass parameters in the matrix M?r are unbounded. If they are very large 
compared to the electro-weak symmetry breaking scale, v, the right-handed neutrinos de- 
couple from the light spectrum. They can be integrated out giving mass to the light 
neutrinos through the effective dimension-5 operator ||19||, 



5 = It [Y u M^Yj] a ^ <£ l{. (2.4) 

which induces, upon spontaneous symmetry breaking, a Majorana mass to the SM neutri- 
nos of the form 



ra v = v 



[Y U M^Y?] . (2.5) 

This is the famous seesaw mechanism j| §, §, f|, that can lead to small neutrino masses, 
even if their Yukawa couplings are of 0(1), provided there is a large hierarchy between v 
and Mat. 
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If the Yukawa coupling of i>l and vr fermions with the Higgs field are, however, of the 
size of the electron Yukawa coupling or smaller, Mn = O(l) TeV or less and new interesting 
signatures are possible. If Mat = 0(100) GeV, direct production of right-handed neutrinos 
through Higgs decay can take place at the LHC fEjj, |2|, |23|. For M N G [0.1, 10] GeV, 
right-handed Major ana neutrinos induce lepton- violating processes in meson decays (for a 
recent study see ^4|). In the keV range right-handed neutrinos can be interesting dark 
matter candidates (for a recent revival see p5|, 26]) or in the eV range, they could contribute 



to explain the LSND-MiniBooNE anomaly [12, [L4|] and affect cosmological observables, such 



as the abundance of light elements, the anisotropies in the cosmic microwave background 



or the large scale structure (for recent analyses see [27 



In this paper we concentrate only on the constraints imposed by oscillation data which 
can be relevant for Mn not much larger than a few eV. 



2.1 Physical parameters 

The number of physical parameters can be determined in different ways. In the following 
we use the method proposed in pOfl . If R/I is the number of real/imaginary parameters in 
the generic Yukawa matrices, R s /I s is the number of real/imaginary parameters defining 
the symmetry group elements of the lepton sector when the Yukawa matrices are switched 
off, and R r jl r is the number of parameters describing the group elements that survive 
for non-zero Yukawa matrices, the number of real/imaginary physical parameters should 
satisfy the following relation: 

Rphys = R — {R s — R r ), 

Iphys = I — {I s — I T ). (2.6) 

We have two separate equations for the real and imaginary parameters: the first ones 
become either masses or angles and the latter become phases. 

Let us do the counting in general for families of active and ur families of steriles. 
The charged lepton Yukawa matrix is ni x n^, while the charged- neutral one is rectangular 
n L x nR. The mass matrix of the Major ana neutrinos is a symmetric square matrix ur x hr . 
The total number of real and imaginary parameters is, therefore: 

R = I = nL { nL + nR ) +nR ^±A. (2.7) 

If all those matrices would vanish, the symmetry group of the lepton sector would be 
U(ul) for the left doublets, U{ur) for the Majorana neutrinos and U(ji£) for the right- 
handed charged fermions. Therefore the number of real and imaginary parameters of the 
corresponding group elements is: 

P = n L {n L + l) + n ^ n l + l \ (2.8) 
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There is no remaining symmetry when the masses are non-zero R r ,I r = 0, and therefore 
the number of physical parameters is: 



Rphys = n L (n R + 1) + n R = n L n R + n R + n L 
Iphys = n L (n R - 1). 



(2.9) 



Since the number of zero modes in the mass matrix is generically til — n R , the number 
of mass parameters is ul + n R + n R = ul + 2n R for > n R (which should be real 
parameters). Therefore: 



It is important to stress here again the difference between the model we are considering 
and the phenomenological 3 + N s models. For example the popular 3 + 1 model has six 
physical angles and four mass eigenstates. In order to have that number of physical angles 
in our case we need to have at least n R = 3, but that choice would generically imply six 
mass eigenstates. In order to recover only four we would need to have some degeneracies, 
that is, some Weyl fermions must pair up into Dirac neutrinos, which typically reduces also 
the number of physical angles. For example, we could get the 3 + 1 model from n R = 5 and 
imposing a global lepton number symmetry (there are probably other choices to obtain the 
same result). 

A more systematic way to constrain models with singlet fermions is to classify them in 
increasing order of complexity, according to the number of extra fermionic field degrees of 
freedom and the global symmetries, in contrast with the classification based on the physical 
spectrum that might represent different choices of the former type. In this work we will 
concentrate on the two simplest cases: n R = 1 and 2. The number of physical parameters 
and the generic spectrum for the different choices of global symmetries are summarized in 



In the absence of global lepton number symmetries, the case n R = 1 gives a spectrum 
of two massive and two massless fermions. There are two physical angles and no CP 
violation. If a lepton number symmetry is imposed, there is only one charge assignment 
for the extra field that allows a renormalizable coupling to the SM neutrino fields (L = 1) 
. In this case, the two massive Majorana fermions are degenerate and form a massive 
Dirac fermion, while two Weyl fermions remain massless. Clearly the model with n R = 1 
with a global lepton number symmetry is ruled out, as it cannot accommodate the solar 
and atmospheric oscillations. On the other hand the model without the global symmetry 
contains in principle sufficient parameters (two mass eigenstates and two angles) to explain 
both oscillations lengths. 

The model with n R = 2 and no lepton number symmetries, gives rise to a spectrum 
including four massive and one massless neutrino. There are also four physical angles and 
three CP violating phases. 

x Any charge assignment that forbids couplings of the singlet fermions to the SM fields is of course 
uninteresting. 




n R (riL - 1) n L > n R 




(2.10) 



Table 2.1. 
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Table 1: Spectrum and number of independent angles and phases for the models with nji = 1,2 
without and with global lepton number symmetries. The second column shows the lepton number, 
L, charge assignments of the extra singlets., L; Only charge assignments were none of the extra 
singlets gets completely decoupled are considered. 



Simplifications also occur when lepton number symmetries are imposed. For ur = 2 
there are two choices for the lepton number charge assignments that allow renormalizable 
couplings between the extra singlet fermions and the SM neutrinos. One obvious choice 
is to give both of the sterile fields lepton number charge +1. In this case, the spectrum 
degenerates into a massless neutrino and two massive Dirac neutrinos. The number of 
physical angles gets reduced to three and there is only one physical CP phase. Obviously 
this choice is as good as the standard three-neutrino mixing model to accommodate existing 
oscillation data. 

The other choice for the charge assignments is to give charge +1 only to one of the extra 
fields and -1 to the other. In this case, the spectrum consist of three massless neutrinos and 
one massive Dirac one. The number of physical angles is reduced to two and there is no 
CP violation. This model with just one mass cannot explain oscillation data. However, a 
small perturbation that breaks the lepton number symmetry is again as rich as the generic 
case of tir = 2, but with some strong hierarchies, naturally preserved by the approximate 
lepton number symmetry. This is the minimal flavour violating seesaw model considered 



in |3lJ (see also |3^]). For ur = 3, there are many more possibilities, listed in Table 2T 
3. Parametrization 

Consider the generic mass matrix corresponding to the model with three left-handed neu- 
trinos and ur right-handed Majorana neutrinos: 



M v = vt« f (3-1) 




where Mn is a tlr x tlr matrix and Y v is a 3 x ur matrix. On the other hand, without 
loss of generality, the Mn matrix can be chosen as a diagonal matrix. For simplicity we 
will restrict to the case where all entries are real, so that there is no CP violation. 
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The first step will be to do a rotation to bring the block Y v vj\pi to a minimal form 
where the rows 1, .., 3— ur can be set to zero and the remaining squared block is off-diagonal 



Mj = n T M^ 



m 
m T M' N 



where 




U T Y U ^=W = m, M' N = W T M N W 
v2 



(3.2) 



(3.3) 



and 



• n R = 1 



m = 




m D 



V2\ 



Y,y v I, U = U 23 (6 23 )U 13 (e 13 ), W = l. 



(3.4) 



a=l 



• n R = 2 



m = 



( \ 

m D - 
\m D + J 



i/ = *7 23 (#23) u 13 (e 13 ) u 12 (e 12 ), w(e 45 



^ ^al — ^a2^j + ^ ^al^a2^ 



1/2 



(3.5) 
(3.6) 



^23; C7i3, C/12 are 3x3 orthogonal matrices corresponding to rotations around the 1, 2 and 
3 axes and the 2x2 W matrix depends only on one extra angle, #45. This procedure can 
be easily extended to the case ur = 3, in which case W would depend on three angles. 

Note that the number of real free parameters of M! v in the new basis (the angles 
included in the matrices U and W and the mass parameters in m and Mat) complete 
the expected total number and provide therefore a complete parametrization. For the CP 
phases, it can be shown that one of them can be absorbed in a <5-type phase in U, and 2/5 
additional phases can be absorbed in W for hr = 2/3 respectively. The latter disappear 
in the Dirac case, Mjy = 0. 

We have not yet brought the full matrix to a diagonal form, but the necessary additional 
rotation angles must be calculable in terms of the elements of M.' v . We first note that the 
Dirac case, Mn = 0, is trivial for any ur, because the final diagonalization is just a trivial 
rotation of angle 7r/4 (i.e. Majorana to Dirac basis) in the subsectors 34 for ur = 1, in the 
sectors 34 and 25 for ur = 2 and 34, 25, 16 for hr > 3. 

Let us consider the diagonalization in the remaining cases. 

3.1 3+1 Model 

For the simplest case ur = 1, the final diagonalization requires a 2x2 rotation on the 34 
subspace: 

M v " = n' T M„'n' = diag(0, 0, A_ , A+), (3.7) 



-7- 



where 



and 



A_ = i (m - ^M 2 + Am 2 D ^ , \ + = 1(m+ ^M 2 + Am 2 ^ (3.8) 



(\ \ 

10 

cos 6*34 sin #34 
\^0 — sin 6*34 cos #34 J 



, sin 6*34 = 



M 



2^M 2 + 4m 2 D 



1/2 



(3.9) 



The full diagonalization matrix is therefore Q f2', where the independent parameters are 
chosen to be (tod, M, #13, #23). 

Obviously we could also use the completely general phenomenological parametrization 
with a 4 x 4 mixing matrix, which relates the mass and flavour eigenstates, U m - lx . Such 
matrix is given, in terms of our parametrization by 



fAmix — ^ ^ • 



(3.10) 



3.2 3+2 degenerate case 

If the Majorana mass matrix is proportional to the identity, Mjy = MI, the matrix M' N = 
Mjsf. The diagonalization of the full matrix involves two independent 2x2 diagonalizations. 
We easily get 



Mj' = diag(A , Ai, A 2 , A 4 , A 3 



(3.11) 



with 



A 
Ai 



A 2 = ^-W(^ 



= 




M 




~ T 




M 




~ T 




M 






+ 


~ T 


M 






+ 


~ T 



MY 2 



to 



D+ > 



M 



+ to 



D- ' 



M 



m 



D+ > 



(3.12) 



where f2' is a rotation on the sector 34 and 25 with angles 

M 



sin #34 



sin t/ 2 5 



2y / M 2 + 4m^_. 



1/2 



M 



2,/M 2 +4to 2 > 



1/2 



(3.13) 
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The free parameters in this case are (M, m^- ,m£>+, #12, #13, #23)- Notice that here the ro- 
tation through the matrix W is unphysical and consequently #45 is not a physical parameter 
in the degenerate case. 

It is easy to see how the Dirac and seesaw limits are obtained and how the parametriza- 
tion reduces to the standard one in the three-neutrino mixing scenario. In the Dirac limit, 
the square neutrino mass matrix is Y y Yjv 2 /2 = Umm T U T = f/Diag(0, m 2 D _ , m 2 D+ )U T 
and therefore the matrix U is just the PMNS matrix. In the seesaw limit, the same is true, 
since the light neutrino mass matrix is Y v M^Yjv 2 j2 = M _1 f7Diag(0, m 2 D _ , m 2 D+ )U T for 
M N = MI. 

The completely general phenomenological parametrization with a 5 x 5 mixing matrix, 
where the relevant matrix elements will be U a i, with a = e, n, r, .. and i = 1,2, 3, 4, 5 would 
be, again, 

u mix = n n'. (3.14) 

3.3 3+2 general case 

The final diagonalization of M! v involves solving the eigensystem of a 2n^ x 2n^ matrix, 
which is not possible analytically in the general case for hr > 1. 

A simplification takes place in the seesaw limit (Y u v Mjv), but in order for it to be 
simple it is necessary to first block diagonalize Ai u to separate the light and heavy sectors. 
We can do this by taking the orthogonal rotation 




where the mixing 9 is small in the seesaw limit. At leading-order 
which satisfies 



Y U ^=M^\ (3.16) 



K=(V'V /r^fV) (3-17) 

" \e T i ) \yJv/V2 m n )\-0 t i) 1 ; 

f-Y v M N l Yjv 2 /2 \ 

" V M N (i + o{e 2 ))) ■ 

(3.18) 

We can finally perform a rotation to diagonalize the light sector, which can be done by 

M'l = n' T M' u n' n '=(o/)' (3 - 19) 

where U can be chosen such that 

U T Y u M N 1/2 Wv/V2 = m 1/2 . (3.20) 
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U can be identified with the usual PMNS matrix (up to phases) while W is a 2x2 orthogonal 
matrix that depends on an extra mixing angle that we can call #45. m 1 / 2 contains the two 
non-zero light neutrino masses: 



m 1 / 2 



( \ 



(3.21) 



We can choose as mass parameters the two masses in m 1//2 and the tir heavy masses in 
M/v- In addition, we have the 3 standard angles in £/(#23> #13 j #12) and the extra angle of 
W{9^). It is important to stress that for the CP conserving case eq. ( 3.2CQ is the usual 



— 1/2 

generic bidiagonalization of the matrix Y U M N while in the general case, it should be 
understood as a definition of the orthogonal matrix W . This parametrization turns out 



to be equivalent to that of Casas-Ibarra [11], where the matrix W is what they call the 
matrix R. 

Note that the matrix f2 can now be reconstructed from these parameters since 

6 = Um 1/2 W T M N 1/2 . (3.22) 
The phenomenological parametrization in this case is 

U mix = nn'~( U nT ~ T 9 ) + O(0 2 ), (3.23) 




and therefore {U m \ x ) a i = a i for a = e, [a,t and i = 4,5. 

Obviously this parametrization is not the same as the one discussed before. However, 
note that in the degenerate case the matrices U and U coincide. This parametrization 
has the advantage to be more physical in the seesaw limit, but cannot be extrapolated 
to the quasi-Dirac limit. For a full exploration of the parameter space that works in the 
quasi-Dirac and seesaw limits, the general parametrization of eqs. (|3.3j ) is more appropriate. 

4. Constraints from oscillations on 3+1 model 

We start by considering the 3 + 1 model. Although the model is very constrained, it has in 
principle sufficient parameters to fit two mass splittings and two mixing angles. Therefore, 
it is interesting to understand to what extent the model can fit oscillation data. 

We have explored the full parameter space trying to fit the data from long-baseline 
accelerator (MINOS ||) and reactor experiments (KAMLAND @, and CHOOZ||), 
since they provide the most precise information on the oscillation frequencies. We carry 
out the analysis of these experiments as explained in [36, p7], and we refer to those papers 
for details. 

On the left plot of Fig. |], we show the contours Xexp = x|m + ^-21 on the plane of 
the two mass parameters (mu, M) after minimizing in the two angles. Xsm corresponds to 
the best fit of the standard three-neutrino oscillation scenario, while Xexp corresponds to 
the minimum in the fits of the ur = 1 model to the experiments (exp): KAMLAND (red), 
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MINOS (green) and CHOOZ (grey). As expected MINOS and KAMLAND allow a very 
constrained region of parameter space that deviates very little from the lines where one of 
the mass splittings equals Am„ tTO and Am^ ( respectively. The two lines intersect only in 
the two points indicated by circles, for the normal (NH) and inverted (IH) hierarchy. On 
those points, the atmospheric and solar frequencies can be accommodated simultaneously, 
but the values of the angles that are required to fit the data are nevertheless incompatible, 
as shown on the two right panels. 

The conclusion is therefore that the 3 + 1 is excluded. Note however that the com- 
bination of the three experiments is essential. This is in agreement with the qualitative 
analysis of ref. fl38j . 




Figure 1: Left: Contours corresponding to xixp = X 2 \cxp +9.21 from the fits of the 3+1 model to 
KAMLAND (Red), MINOS (Green) and CHOOZ (grey) on the plane (m D ,M) (after minimizing 
in the angles). The solid lines correspond to the values of tod and M such that one of the mass 
square differences in the spectrum coincides with the solar or atmospheric splitting. The circles 
indicate the values of rriD,M where the solar and atmospheric splittings coexist for the normal 
(NH) and inverted hierarchy (IH) . Right: same contours on the plane (#13, #23), for the values of 
tod, M fixed at the intersection points NH and IH. 

5. Constraints from oscillations on 3+2 model 

For the case tlr = 2, the number of parameters is quite large to do a full-fledged fit as has 
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been done in the 3 + 1 model. We will consider the simpler degenerate limit, where there 
are only three independent masses and three angles and furthermore consider the CP con- 
serving case. In order to identify the allowed regions we will follow the strategy of looking 
for solutions where at least one of the mass splittings corresponds to the atmospheric one 
and another one to the solar. Given that both splittings are well established by data, this 
method should be rather robust in pinning down local minima. 

More concretely we consider all possible combinations ijkl such that 

j A? — | = Am^ tm = 2.5 x 10" 3 eV 2 , |Aj£ - \f | =Am s 2 ol = 8x 10 _5 eV 2 , i,j,k,l = 0, ..,4 

(5.1) 

These two equations imply relations between M, m 2 D + and tn 2 D - that we solve numerically 
for ra^y and m~jy as a function of M. In principle there could be up to 90 combinations, 
but we restrict to those that correspond to rn? D+ > fn 2 ~ ) _, since they have been ordered 
in this way. In the quasi-Dirac limit M <C rri£,+ ,m.£)- , most combinations are possible, 
while in the opposite limit M 3> m D ± two states decouple and therefore only the standard 
choices of the 3v scenario survive (two for the normal hierarchy and two for the inverse 
one). We can therefore classify all the possible solutions by the hierarchy they represent 
in the quasi-Dirac limit. We find the following, qualitatively distinct, spectral patterns 
depicted in Fig. ||. 

• Type I: Quasi-Dirac Normal Hierarchy 

For small M, we could have the solar splitting between eigenstates and 1 or and 
3, and it is therefore related to m D - ~ J Am 2 ol . The atmospheric splitting could be 

either between and the heavier states 2 or 4, ie. m D + ~ y Am 2 tm . But it could 
also be between the states 1/3 and 2/4. There are in total 12 combinations that 
differ only by small perturbations of 0(Am 2 ol /Am 2 tm ). Typically these solutions 
have other mass splittings that are much smaller than both solar and atmospheric 
for sufficiently small M. As M increases, most of these solutions disappear except 
the two that involve only the states 0,1,2. 

• Type II: Quasi-Dirac Inverse Hierarchy 

This case is like type I but representing an inverted hierarchy. This happens when 
the atmospheric mass splitting occurs between states and 1/3 or 2/4, while the 
solar one does between states 1/3 and 2/4. In this case there are 16 combinations 

in total, which have m D ~ ~ \J Am 2 tm and rn 2 D + — — \J ^ m s i- Again most of 

these solutions do not survive for large M, only two do. All the additional splittings 
are not larger than those measured. 

• Type III: Quasi-Dirac Degenerate 

In this case, the solar splitting is the one between states 1 and 3, while the atmospheric 
one is the splitting between states 2 and 4. This solution is only possible for not too 
large M and it is rather peculiar because it has other splittings that are larger than 
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Figure 2: Spectral structures that contain the solar and atmospheric splitting in the quasi-Dirac 
region. 



those measured, since typically m D +,m D - S> a/ Am^ tm , y Am^. It is extremely 
interesting if it gives a good fit to the data. 

• Type IV: Quasi-Dirac Solar Degenerate 

This is a mixed situation between the quasi-Dirac and quasi-degenerate cases. There 
are two possibilities. In type IV-a, the solar splitting is the one between states 1 and 
3, while the atmospheric one is driven by mrj+, that is the splitting between states 
and 2 or 4 or alternatively 1 or 3 and 2 or 4. There are 6 possibilities of this type. 
Type IV-b corresponds to a solar splitting associated to the states 2 and 4. Then the 
atmospheric must be related to m D -, so it is the difference between states and 1 
or 3. Two combinations exist of this type. 

• Type V: Quasi-Dirac Atmospheric Degenerate 

In this case, the atmospheric splitting is the one between states 2 and 4, while the 
solar one is between states and 1 or 3. There are two possible combinations of this 
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type, which correspond to a situation with mp- ~ JAmL and m£)+ 3> 771d- 



There are some more solutions that only occur in the intermediate region of M where no 
clear hierarchies can be established. These however give a very poor fit to the data, so we 
will not discuss them further. 

5.1 Fits to reactor and accelerator experiments 

We started by considering the fits to long-baseline reactor and accelerator experiments 
(LBL), that is including Chooz, KamLAND, MINOS appearance and disappearance for 
neutrinos and antineutrinos, that in principle have the most precise information about the 
two mass splittings. We find that solutions of all the types mentioned provide a good fit 
to the data, as good as the standard three neutrino scenario. That this holds is obvious in 
two limits: the Dirac limit M — > 0, where the spectrum matches a three-neutrino mixing 
scenario with one massless neutrino, and, for those solutions that survive, the seesaw limit 
M 3> m D ± , where the spectrum also corresponds to a three-neutrino mixing scenario with 
one massless neutrino, since the two heavier states decouple. The survival of other solutions 
with M between the two limits is however non-trivial since they imply the presence of other 
mass splittings that could affect oscillations. 

In Figs. || we show the minimum x 2 after minimization in all the angles, as a function 
of M. We only show those solutions that give a good fit for some value of M, which can 
be classified according to the hierarchies above. 

The left upper plot corresponds to the solutions of type I . As expected the x 2 increases 
significantly with respect to the Dirac- limit value for intermediate values of M. From the 
point of view of LBL data we see that values in the range 5 x 10~ 4 eV < M < leV give 
poor fits to the data. There are however also some values of M in this range that give 
acceptable fits to the data for some of the solutions, that need further scrutiny. 

A slightly different situation is found for the solutions of type II, which for small/large 
enough M correspond to an inverse hierarchy. Here LBL data can exclude a slightly wider 
region 10 -4 eV < M < leV and we do not see any local minimum that provides a good fit 
within this range. 

The quasi-degenerate and mixed solutions are shown in the lower plots. In this case 
the solutions only exist for sufficiently small values of M < 5 x 10 -3 eV, as shown. The plot 
on the left corresponds to types III and V, while the one on the right corresponds to types 
IV. The type III gives a rather good fit to the LBL data in all the region, while solutions 
of type V only for M < 5 x 10 _4 eV or for M ~ 0.005 eV. Some of the solutions of type IV 
instead give a good fit in most of the range where they exist: 10 _5 eV < M < 10 _3 eV. 

In Figs. |] we show the minimum Xlbl on ^ ne pl anes defined by each pair of mixing 
angles for the solutions of type I and II as examples. Typically one minimum is found, but 
there are octant degeneracies, as expected. The solutions can be easily identified by the 
very different positions of their minima on these planes. 

Note that we consider the same ordering of states in all cases and therefore the so- 
lar and atmospheric mass splittings are labelled differently for the different cases. As a 
result, only in the normal hierarchy cases the values of the angles 612, 613 and 623 are in 
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Figure 3: Minimum Xlbl fr° m the fits to KamLAND, MINOS and Chooz as a function of M(eV). 
The solid(dashed) horizontal line corresponds to Xsm(Xsm + 9) of the standard Sv scenario. Up 
left: solutions of type I. Up right: solutions of type II. Down left: solutions of types III and V. 
Down right: solutions of type IV. 

one-to-one correspondance with those in the standard parametrization of the three neu- 
trino scenario. In the inverse hierarchy case, the angles are quite different, although the 
oscillation probabilities in vacuum will be equivalent. For example, the normal hierarchy 
for the angles (#12 = S oi,0i3 = 0,023 = a tm) is equivalent for the inverse hierarchy to 
(#12 = ttA 013 = 0soi, 023 = 7T - a tm), as can be seen in Figs. |. 

In summary, LBL data are able to exclude most solutions in an intermediate region 
of M, but some peculiar solutions do survive. As expected, the quasi-Dirac type I and II 
solutions survive for all M < Mmax, and those involving only the lighter 012 states also 
survive in the seesaw limit, M > M^f n . 

In the following we will further constrain the intermediate solutions using solar and 
atmospheric data, and also refine the determination of the quasi-Dirac and seesaw bound- 
aries: Mmax and Af£f n . 

5.2 Solutions for intermediate M 

Solar data give additional constraints on solutions in the intermediate region (we refer to 



1 36] for details of the data analysis). The crucial piece of information provided by solar 



data is of course the neutral current measurement of SNO [39]. In Fig. [| we compare the 
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Figure 4: Xlbl contours (l,2,3cr) on the planes (#12, #13), (#12, $23) and (#13,023) at fixed M as 
shown, and for solutions of type I (up) and II (down). 



X 2 of the LBL fit and that of solar data to those of the standard 3u scenario (horizontal 
lines) for M in the meV range, 10 -3 -10 -2 eV. All solutions that survive LBL are essentially 
excluded by solar data. The same is shown in Fig. || for the ceV region, 10 -2 -10 -1 eV. In 
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Figure 5: Xlbl (green triangles) from the fits to KamLAND, MINOS and Chooz and Xsun ( rec ^ 
diamonds) as a function of M(eV). The horizontal lines corresponds to the Xmin °f the standard 
3v scenario for the LBL (green) and solar (red) fits. The dashed lines correspond to a shift by 
A X 2 = 9. 
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Figure 6: Xlbl (green triangles) from the fits to KamLAND, MINOS and Chooz and Xsun ( re d 
diamonds) as a function of M(eV). The horizontal lines corresponds to the Xmin °f the standard 
Zv scenario for the LBL (green) and solar (red) fits. The dashed lines correspond to a shift by 
A X 2 - 9. 



this small region around M > 10 2 eV survives both fits for types I and V. For 

those cases, we have considered the fit to atmospheric data [40] and the result is that none 
provide a decent fit to atmospheric data and are therefore excluded. 



5.3 Limit on quasi-Dirac neutrinos 

Solutions of types I and II provide a perfect fit to the data for M < Mmax, which bounds 
the quasi-Dirac region. Including solar data is essential in constraining Mmax to much 
lower values that those obtained from just fitting LBL data, for which Mmax ~ 5 x 10 -4 
eV (see Fig. |3|). The upper bound Mmax is therefore constrained exclusively by solar data. 
The result is shown in Fig. [7| for the solutions that have a Dirac limit, i.e. those of types I 
and II. We obtain a value of Mmax — 10 -9 , 10 -11 eV for the NH, IH solutions respectively. 
In the last section we estimate analytically such extremely small limits. 



5.4 Mini-seesaw region 

For M sufficiently large, the states 3 and 4 decouple from the spectrum and only four 
solutions provide good fits to LBL data (two corresponding to NH and two to IH). The 
value at which this seesaw limit sets in is M^ n and it is essentially determined by LBL 
data, as shown in Fig. g. We find M^f n ~ 0.6 eV, 1.4 eV for the NH, IH solutions, 
respectively. 

Neither solar nor atmospheric data add any further constraint on M^^. The mini- 
seesaw solutions provide a perfect fit to solar and atmospheric data for larger values of M. 
Near the seesaw threshold the mini-seesaw solutions that agree with solar, atmospheric 
and LBL data have a hierarchy of the form shown in Figs. p|. 
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Figure 7: Minimum X%un fr oin solar data as a function of M(eV) in the quasi-Dirac region for 
solutions of type I (NH) and of type II (IH) . The solid horizontal line corresponds to the minimum 
X 2 and the dashed one to a shift by A% 2 = 9. 




Figure 8: Minimum Xlbl as a function of M(eV) in the mini-seesaw region for solutions of type I 
(NH) and of type II (IH). The solid horizontal line corresponds to the minimum \ 2 and the dashed 
one to a shift by A\ 2 = 9. 



6. Other bounds in the mini-seesaw region 



6.1 LSND/MiniBoone and reactor 

Obviously, an interesting question remains concerning the neutrino anomalies: the old 
LSND result recently reinforced by the anti-neutrino appearance data in MiniBoone 
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Figure 9: Spectra of neutrino mass eigenstates in the mini-seesaw region for M ~ M^ n , M 
0.6eV for the normal (left) and M ~ 1.4 eV for the inverse hierarchy (right). 



1 1*7 , 18], the low-energy excess of MiniBoone and the recently discovered reactor neutrino 



anomaly pi] , E2fl . A possible explanation of such anomalies is the existence of extra sterile 



neutrinos in the eV range [13, ||, [14|]. It is then an obvious question whether our 3+2 model 
for M > M^ n can account for any of those anomalies. 

In [|l4|] a recent reanalysis of the 3+2 phenomenological scheme shows that there is a 
significant improvement in the goodness-of-fit with respect to the analysis previous to the 
reevaluation of the reactor fluxes, and to the publication of the MiniBoone antineutrino 
data. Nevertheless, there is still a tension between the LSND/MiniBoone appearance and 
the muon disappearance experiments. 



According to flL4| , a good fit to the data is obtained with: 

1. Two neutrino mass differences around 1 eV (more precisely around 0.5 eV and 0.9 
eV). 

2. In the phenomenological parametrization, large enough |£7 e 4|, |£/ e5 |, {U^l and \U IJ ^\, 
around 10 _1 . 

3. CP-violation. 

Since the parametrization for the phenomenological 3+2 model contains all the possible 
physical parameters when there are 5 mass eigenstates, our model can be described by the 
same parametrization although there are correlations between the parameters, in partic- 
ular between the masses and angles. In particular for the degenerate case, according to 
eq. ( p. 23 ), we have 



(C^mix)e4 = S13S34, 

(C^mix)e5 = C13S12S25, 

( fornix )/i4 = C13S23S34, 

(^mix)jtt5 = (C12C23 - Si 2 Si 3 S23)s 2 5, (6.1) 
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with 



S25 ~ rn D - /M, s 34 « -m D +/M, (6.2) 

for M in the eV range. 

Concerning the massive states, their masses are: (~ m 2 - ) _/M, ~ m 2 D+ /M,^ M, ~ M). 
Therefore, after taking into account the results from the above sections, the only parameter 
we can play with is the Majorana mass M, and therefore in the degenerate case it is not 
easy to accommodate two distinct eV masses. On the other hand the matrix elements 
U e 4/ 5 are intriguingly in the right ballpark. 

In order to have a more explicit prediction we have evaluated the P ee and P e/Ji oscillation 
probabilities in the mini-seesaw regime analytically via a perturbative expansion in the 
small parameters 

( m D-\ 2 ( m D+\ 2 f(t „. 

e = hsr) ' e + = hr) ■ (6 - 3) 

At second order in e± we get for the ten mass differences: 

Am 2 ,- = 0(e 2 M 2 ), Am| ~ Am| ~ M 2 (l + 0{e 2 )), Am^ 4 ~0(eM 2 ), i, j = 1, 3.(6.4) 
At the same order in e± the mixing angles are: 

sin 2 #25 = e_ — 3e 2 , sin 2 #34 = e+ — 3e+, 

and finally the oscillation probability is 

'M 2 V 



P e f = 4sin 2 ( — \ (e 2 + A u + 2e+e.A 54 + e 2 _A 55 ) + 0(4), (6.5) 



where 



J**™ = 4sin 2 ( ^ ) ^lAnh + 0(4), (6.7) 



^44 = cos 2 #13 sin 2 #13 sin 2 #23, 

^54 = cos 2 #13 sin #12 sin #13 sin #23 (cos #12 cos #23 — sin #12 sin #13 sin #23 ) , (6.6) 
^55 = cos 2 #13 sin 2 #12 (cos #12 cos #23 - sin #12 sin #13 sin #23) 2 . 

In the normal hierarchy case, #13 is small and will introduce a suppression of the ^44 and 
^54 coefficients with respect to For vanishing #13 we get: 

'M 2 V 

with Anh = c i2 c 23 s i2 = 1 sm2 2#soi cos 2 # a tm ~ 0.1. On the other hand, in the inverted 
hierarchy case, #12 ~ vr/2. In the limit #12 = vr/2 we get: 

P e f^ = 4sin 2 (^) (e+ - e „) 2 ^ + 0(4), (6.8) 

with Ajh = c 2 3S 2 3 s 2 3 = A sin 2 2# so i sin 2 # atm ~ 0.1. The inverted hierarchy gets therefore 
an additional suppression with respect to the normal hierarchy case due to the partial 
cancellation between e + and e_. It is easy to see that 



{( - l -'- l ^]^ L ^- H - (6-9) 
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and, therefore, a suppression of the order of the solar to atmospheric mass ratio is to be 
expected in the inverted hierarchy case with respect to the normal hierarchy one. Unfortu- 
nately this accidental suppression for the IH is responsible for the fact that the leading-order 
result of eq. Q6.8| ) is not sufficiently precise and higher orders are relevant, but in any case 
the effect is too small. The numerical expressions for P e/i for L = 541 m in the range of 



energies relevant for MiniBooNE are shown in Fig. 10 for the NH and IH, compared with 



the best fit result of [14]. The strong suppression for the IH case is clearly seen. The 



degenerate case cannot accommodate the LSND anomaly. 

A different situation is found in reactors, which actually set stronger constrains on 
M^f n . The reason is that the effects on disappearance are only linear in e and not quadratic. 
In the same perturbative expansion as before we find: 

Pef = 1-4 + e_c? 3S f 2 ] sin 2 

M 2 L 

P%? = 1-4 [e+c 2 13 s 2 23 + e_(ci 2 c 23 - s 12 s 13 s 2 3) 2 ] sin 2 (6.10) 

This is the origin of the well-known tension between the LSND appearance signal and the 
disappearance constraints, in the electron channel by the Bugey-3 experiment [|43|] and in 



the muon channel by the CDHSW experiment [44[. 



In Fig. |ll| we compare the exclusion plot of Bugey-3 and CDHSW with the prediction 
for the effective mixing angle and frequency in the disappearance channels in the degenerate 
case (note that the two are related). We can see that there is a very small additional 
exclusion of the allowed region for the IH, which sets 

M*f n ~ 1.6 eV, (6.11) 

while there is no further constrain for the NH case. Note that these bounds would be 
relaxed if the new computation of the reactor fluxes is taken into account. 

6.2 Other bounds 

It is clear that important constraints on the mini-seesaw region could also come from 
tritium /3-decay [E3, 46], which is sensitive to the combination 2 



/^|C/ ei | 2 m 2 i <2eV, (6.12) 



and from cosmology. The amplitude of the neutrinoless double /3-decay process on the 
other hand approximately vanishes for M <C 0(100) MeV and therefore cannot constrain 
further the mini-seesaw region (see pSfl). 

Naively one would expect that such bounds will be in the form of a maximal value of 
M^f x in the mini-seesaw region. 



2 It is assumed that all m, <C AE, where AE is the energy resolution near the end-point. When non 
degenerate neutrinos are considered the analysis is more complicated |l7| and cannot be cast as an upper 
bound on the combination m e . 
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Figure 10: The P^ e antineutrino oscillation probability at L — 541 m as a function of E v for the 
NH (red, M = 0.6 eV) and IH (blue, M = 1.4 eV). The best fit result of @ is labelled KMS. The 
dashed lines correspond to the perturbative results. 




Figure 11: Exclusion region at 99%CL from the experiments Bugey-3 (left) and CDHSW (right). 
The lines correspond to the expectation as a function of M in the degenerate case. The intersection 
of the IH curve and the Bugey-3 exclusion region is at M ~ 1.6 eV. 

Concerning tritium /3-decay, at leading order in the e expansion we get 



m e - V S?3 m D+ + C 13' S 12"4>- ■ ( 6 - 13 ) 

For the NH we therefore have 

m e ~ yj S 2 13 ^/Aml tm + Am^^M ~ 0.05 ^M{eV)^- (6-14) 
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For the IH instead (#12 ~ 7r/2 and #13 ~ so i): 



m e ^y [s 2 13 + c 2 3 s 2 2 )^Am 2 atm VM ~ 0.22 v / M(eV)eV. (6.15) 

Therefore, the present bound for the degenerate case gives M^f x < 82 eV (IH) and M^f x < 
1.6 keV (NH). These estimates are however too naive ||, because for such large values of 
M, the heavy states are far from the end-point and the sensitivity to the heavy states 
cannot be cast into the form of their contribution to m e in eq. ( ljj) |7fl . In fact, a more 



carefull analysis along the lines of || shows that there is no bound on M from the present 
measurements of tritium experiments, neither for the NH not for the IH. 

More stringent bounds are expected from cosmology, that imposes stringent constraints 
on the number of thermalized relativistic species at the time of nucleosynthesis (BBN), and 
at recombination (CMB) p9[ . Furthermore the large scale structure is already sensitive to 
neutrino masses in the eV range. Although a detailed calculation is required, both sets of 
measurements will very likely add new information on M^f x maybe closing the eV window 
for this model. 

7. Towards the non-degenerate case 

Although a systematic exploration of the parameter space in the non-degenerate case will be 
studied elsewhere, in this section we want to point out the usefulness of the parametrization 
introduced in sec. ||, which makes it easy to study analytically in the mini-seesaw region 
the oscillation probabilities in the same perturbative expansion introduced before. We 
consider as an illustration the and P^ e probabilities to show that the re# = 2 model 
in the non-degenerate case could explain the LSND/MB/reactor anomaly similarly to the 
solution found in [14]. 

We consider the expansion up to C(e 2 ), where e is any of the ratios of the light to 
heavy mass eigenstates 

e ij =m l /M j , (7.1) 
with rrii the entries in eq. (3.21). These parameters reduce to e± in the degenerate limit. 



The results at leading non-trivial order are 

P e f = 4|0 e4 | 2 |M 2 sm 2 ^ +4|tf e5 |X 5 | 2 sin 2 ^ 



(M$-M?)L \ . MfL . M\L 



+ 8|0 e4 M 40e50/,5| COS I ^ \ E + J Sm sin + - . ( 7 - 2 ) 

= 1 " 4|M 2 sm 2 ^ - 4|^ 5 | 2 sin 2 ^ + ... , (7.3) 

where 4> = arg {0e^*^G* e ^^)- The CP-conserving case corresponds to (j) equal tt or 0. We 
introduce the CP non conservation only through (j) as a first step 3 . Using eq. ( [3.22| ) 

e = u(e 12 ,e 13 ,e 2 3)m 1 / 2 w(e 45 fM- 1/2 , (7.4) 



3 A specific parametrization in terms of CP phases in the mixing matrix, will lead generically to phase 
dependence also in the \9 a i\. 
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it is straightforward to obtain 9 a i at the required order. 
The appearance probability is 



= 4 A 44 sin 2 + 4 A 55 sin 2 



r- — — f(M2-M?)L \ M$L m?l 
+ 8 y/A^A^ cos ^ 2 4g + <Aj sin sm 

where, for NH and #13 = 0: 

aNH 22/1/2 1/2 N 

-4 4 4 = £2lSl2'S45 ( e 31 c 45«23 + £ 2 1 Cl2C23«45 



(7.5) 



A NH __ 
^55 — 



2 2 / 1/2 1/2 \ 2 „ N 

e 2 2Si 2 c 45 (^e 22 ci2C2 3 c 45 - e 32 s 23 s 45 J , (7.6) 



and for IH and #12 = vr/2: 

2/1/2 1/2 \ 2 ( 1/2 1/2 \ 2 

^44 = s 23 [ e 31 c 45 s 13 + H\ c 13^45 J (^31 c 13 c 45 ~ e 2 l S 13S45 J , 

/i/H 2/1/2 1/2 \ 2 / 1/2 1/2 \ 2 , v vA 

^55 = S 23 (^22 C 45S13 + € 3 2 C 13«45 J ^22 C 13C45 ~ £32 Sl3«45 ) > ( 7 - 7 ) 

The electron disappearance probabilities are 

P « S |jVff = 1 " 4sin 2 12 ^21^5 ^ ^ + ^4 ^ ^) . ( 7 - 8 ) 

1/2 ,1/2 \ 2 . 2 M 2 L / 1/2 1/2 \2 . 2 M|L" 

e 3 i c 45 si 3 + e 2 [ s 45 ci 3 j sm — — + [e 22 c 45 ci 3 - e 32 s 45 si 3 j sm 



p 55| 
ee \ih 



4E 



It is easy to check that these results coincide in the degenerate limit with those in the 
previous section, as they should. 

Setting the two heavy states to the masses indicated in [jl4|] , it turns out that the angles 
are in the right ballpark for the IH. For the NH on the other hand, they are too small. In 
Fig. [TJ we compare the probability for the best fit of [14]. We find that for 645 ~ 20° and 



4> = 1.627T the results for the appearance probabilities are quite similar and therefore the 
model can also reproduce the LSND/MiniBoone anomaly with a similar level of agreement 
as 1 14]. The signal in disappearance shown in Fig. [l^ is slightly larger so probably there is 
more tension between appearance and disappearance. 

Finally and in the same expansion we obtain for the tritium (3 decay effective mass of 
eq. ( |6.12|) , at the leading order 

(1/2 1/2 \ 2 / 1/2 1/2 \ 2 

m 3 c 45 si 3 + m 2 ' ci 3 si2S 45 J + M 2 \m 2 ' c 13 c i5 si 2 - m 3 ' s 13 s i5 j . (7.9) 



ml 



Fixing the free parameter 845 to zero, a simpler expression which is the same as in the 
degenerate case, eq. ( 6.15Q , is obtained. For the KMS best fit choice of M\ and M2 we 



obtain m P ~ 0.05 eV for NH and m e ~ 0.2 eV for IH. 
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Figure 12: Comparison of the energy dependence of the appearance probabilities P^e for neutrinos 
(right) and antineutrinos (left) at the MiniBoone baseline (L=541m) for the solution KMS [Q (solid 
line) and the non-degenerate case (IH) fixing M\ and M% to the best fit values of [|l4| . The bands 
correspond to varying the free parameter #45. 
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Figure 13: Comparison of the energy dependence of the electron disappearance probabilities at 
the Bugey-3 baseline (L=40 m) for the solution of [O (solid line) and the non-degenerate case 
fixing Mi and Mi to the best fit values of O] . The bands correspond to varying the free parameter 

#45- 



8. Quasi-Dirac neutrinos in the sun 



In this section, we do an analytical study of the flavour transitions of quasi-Dirac neutrinos 
(solutions of Types I and II) in the sun in order to understand the very stringent limits on 
the value of M set by solar data. A similar study has been discussed in [ 50 1 , see also pl|. 



Older studies of pseudo-Dirac neutrinos in the sun were done in |52 



In order to understand the reason for the sensitivity of the sun to tiny splittings in the 
standard three Dirac neutrino scenario, we consider the NH case with #13 = 0, which is a 
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good approximation. The IH works similarly, after the appropriate mapping in the angles 
(in particular, #13 = is equivalent to #12 = 7r/2), so in order to simplify the discussion 
we just consider the NH case in the intermediate results and indicate the differences in the 
final results for the IH. 

Let us first consider the case M = 0, where the standard scenario must be recovered. 
The eigenvalues in matter are: 



(o) 
Mo 



a {0) 
pi 



B 2 



^\A + m 2 D . 



= m D+ 

B 

= — + m 



D+ 



A 2 + 



in 



D- 



2Am 2 D _ cos 26*12 



A 2 + mjj. 



2Am 2 D _ cos 26*12 



(8.1) 



where A = 2y/2G F N e E v and B = 2-^2G F N n E v . Since we are interested in a rough 
estimate, we will neglect N n in front of N e , that is B in front of A, since in the sun this 
is not too bad an approximation. The corresponding eigenvectors, v\°\A), fall in three 
invariant subspaces : the 1, the 4 and 023. The standard 3v result corresponds to the 023 
sector where the usual two- family MSW is found in the sector 02. The MSW resonance 
occurs at Amsw = m \>- cos 2#i2- However, the eigenvalue becomes exactly degenerate 
with /4° } at A = 0. For M = 0, there is an exact level crossing because the levels do not 
mix. When M is non-zero, no matter how small, there is mixing and therefore an additional 
MSW effect in the sector 12 takes place, where the mixing is controlled by M, as we will 
see. 

It is found that only the states vf^ with i = 0, 2 have an electron component, and 
therefore the electron flavour state produced at the center of the sun is a mixture of both 
vf ) \Ao), where Aq corresponds to the value of A close to the center of the sun. In the 
adiabatic approximation, therefore, the probability v e — > v e (ignoring possible effects due 
to Earth matter) is: 



(4Vo))e" 
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We recover the expected Dirac result: 



(8.2) 




12 A » Amsw, 

\ sin 2 2012 Aq < Amsw- 



The result for the IH case is equivalent after exchanging m 2 
The relevant mass scale is therefore in both cases Am s 2 ol . 



d- 



m 



D+~ 



(8.3) 



m 2 D _ and 6*12 ->■ 6*13. 
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Let us now consider a value of M <C mjj- < m,£)+. We can do perturbation theory to 
include the effects of M for any value of A, but for A — > 0, the pairs of states 12 and 34 
are degenerate, and therefore standard perturbation theory breaks down. Non-degenerate 
perturbation theory is necessary here, which amounts to diagonalizing the perturbation 
matrix exactly in the corresponding 2x2 subspaces, treating the rest in standard pertur- 
bation theory. The resulting eigenvalues and eigenvectors are denoted by n\ M ^ an d V \ M ^ • If 
we neglect all effects of 0(M/m D ±), we have that, for large A, the eigenstates go smoothly 
to those at M = 0: 

v\°\A). (8.4) 



lim v^ M \A) 



However, for small A, more concretely for A <C Mmj)± , the M — > limit does not coincide 
with the M = case: 



lim (v[ M) (0j) = - 



M-*0 

lim U M \0)) 



sinfi2 
sin 6*i2 



(8.5) 



while (^ 0) (0)) e = 0. 

In the adiabatic limit, we find 
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.6) 



Note that for energies above the MSW resonance P ee is 1/2 of the Dirac result, so, to the 
extent that the adiabatic approximation is valid, the result for M / is physically very 
different to the Dirac limit. 

A similar analysis in the IH case is a bit more complicated, because in this case the 
problem involves the two pairs of degenerate states and not just one pair as before. The 
final result is 



sin 2 0i3 

2 

o4 



A > A M sw, 



(8.7) 



where the MSW condition is the usual one with the change m 2 D _ — > rn 2 D+ — m 2 -,-. Note 
that in our parametrization, #13 is the solar angle for the IH case. In this case, the result 
is half of the Dirac result both above and below MSW. 

Obviously we expect the limit M — > to be smooth and this means that adiabaticity 
must break down for small enough M. We show that this is indeed the case. 
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8.1 Adiabaticity limit 

Adiabaticity is lost in the propagation inside the sun when 



i/4 M) -^4 M) i 
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where r is the radial distance. For m D ± large compared to the other scales we have 
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(8.9) 



The variation of A can be approximated by 



dA 

— ~ -aA/R & , 
dr 

where a ~ 10 — 15 and Rq is the solar radius. 

The right term in eq. flB.8|) is maximal at the point in the evolution where 
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AM 2 m 2 D . 



At this point the non-adiabaticity condition reads 
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4:\/2R Q m D - ' 

The result for the IH is the same with the change #12 — > #13, but since m^- 



.11) 



.12) 
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for NH and nip- ~ y Am 2 tm for the IH, the value of M for which the adiabaticity limit is 
reached is lower for IH than for NH, by the ratio of solar to atmospheric mass splittings. 
The rough estimates are 



M(eV) < 



10~ 7 x E v (MeV) NH, 
2 x 10~ 8 x E u (MeV) IH. 



U3) 



8.2 Dirac Limit 



Now that we have established that adiabaticity is lost, we must still prove that the Dirac 
limit is reached for such small M. An easy way to see how this happens is to consider the 
sudden approximation. The evolution of the eigenstates shows an abrupt change in the 12 
sector at the condition of eq. 



close to (A), but at this threshold, the eigenstates 1 and 2 change abruptly to be similar 
to those in v\ M \q). The sudden approximation assumes that the physical states do not 
change abruptly and therefore a physical state that is an eigenstate as it approaches the 
transition for A > A t h does not remain in an eigenstate for A < A t h. Oscillations therefore 
occur. 



.11). For larger A the eigenstates evolve smoothly and are 
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At the center of the sun the electron neutrino is a combination of Vq M \Aq) ~ v^ 1 ' (Aq) 
and v^\Aq) ~ v^\Aq) eigenstates. The two very fast decohere and the two components 
evolve adiabatically until the transition where 
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where we have neglected effects of 0(A t h/m 2 D± , M/m D ±). The evolution of the state is 
adiabatic throughout (this is different for the IH case). 

The time evolution of the state after the transition point at io is approximately 
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Using eq. (p 
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the P ee probability in the sudden approximation is found to be: 

Mm, 



sin 2 #i2 cos 2 



ee\NH 



2E V 



-At 



cos 4 6\2 + sin #i2 cos 2 



MSW 



Mm, 



2E V 



-At 



A <A 



MSW- 



To result for the IH is instead 
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When the oscillatory terms are approximately one, the Dirac limit is reached, while the 
regime of fast oscillations reduces to the adiabatic result. In order to recover Dirac we need 
two conditions: 1) non adiabaticity and 2) long enough baseline for unaveraged vacuum 
oscillations with the splitting Mm D -. At the adiabaticity limit (eq. J8.12| )), the oscillation 
length is roughly the solar radius for all energies, and therefore the vacuum oscillations 
are averaged out on Earth. Smaller values of M are necessary to increase the oscillation 
length. These features are precisely found in the exact P ee and P ea shown in Fig. 15 for 
the NH and in Fig. [l5| for the IH. For M ~ 10~ 5 eV, the adiabatic result of eq. (|8.12|) is 
found for all the energies shown. For M ~ 10 _6 eV the vacuum oscillations are seen for 
larger energies while there are averaged oscillations still for the lower energies. For smaller 
values of M, the vacuum oscillations are seen at all energies, and for M < Mmax the Dirac 
limit is obtained. 



9. Conclusions and Outlook 

Probably the simplest explanation of neutrino masses involves the addition of singlet 
fermions to the SM. Such is the case in very different models ranging from Dirac neu- 
trinos, to type I seesaws, inverse seesaw, etc. All these possibilities, that have very dif- 
ferent phenomenological implications for flavour physics, correspond to different numbers 
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Figure 14: For the NH, P ee (red), P ea = P ee + P ell + P eT (blue) at day time as a function of neutrino 
energy for solar neutrinos (B) in the quasi-Dirac region for three values of M = 1CP 9 , 10 -6 , 10 _5 eV 
and normal hierarchy. The dashed curves correspond to the standard 3v solution (Dirac limit) near 
the best fit and the solid lines are the exact results in the quasi-Dirac for the same values of the 
parameters. The production point has been averaged out. 





M=10" 5 eV 



Figure 15: The same as Fig. 14 for the IH. 



of extra species and/or different global symmetries. It is important to study the present 
constraints on models with singlet fermions in increasing order of complexity, where com- 
plexity is measured by the number of extra degrees of freedom (Weyl fermions), tir, and 
not by the physical spectrum, because the latter can depend also on the global symmetries 
imposed. With this perspective in mind, we have considered in the present work the two 
simplest possibilities, that of one or two extra Weyl fermions, ur = 1,2. The first case has 
sufficient free parameters (two masses and two mixing angles) to fit in principle the two 
confirmed oscillations (solar and atmospheric) , but a complete scan of the parameter space 
shows that it cannot fit all the available data from long-baseline reactor and accelerator 
experiments (LBL). The case ur = 2 is equivalent to the standard 3z/- mixing scenario (with 
one massless neutrino) in two limiting cases: the Dirac limit (vanishing Majorana masses 
for the extra fields) and the seesaw limit (large Majorana masses for the extra fields), and 
as such it does provide a good fit to the data in both limits. What happens in between 
these limiting cases is much more complicated. We have thoroughly studied the allowed pa- 
rameter space in this model by requiring that at least one mass splitting corresponds to the 
solar one and another to the atmospheric one, and considered the simplifying assumption 
of degenerate Majorana masses, M. Even though many exotic solutions for intermediate 
values of M can fit very well reactor and accelerator data, they are shown to fail to explain 
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solar or atmospheric neutrino data. We exclude therefore all regions except the quasi- 
Dirac, M < M^ax = KT 9 (1(T 10 ) eV, and the mini-seesaw for M > M^ n = 0.6(1.6) eV 
for the NH(IH) respectively. The upper bound, Mmax, is essentially constrained by solar 
data alone, while the lower bound M^f n is set by LBL and reactor disappearance data. A 
relevant upper bound to the mini-seesaw region could clearly be obtained from cosmology. 

We have discussed in detail the very stringent constraint imposed by solar neutrino 
experiments on quasi-Dirac neutrinos. 

The relevance of the mini-seesaw solutions to explain some of the unsolved neutrino 
anomalies: reactors, LSND, MiniBoone is discussed and discarded in the degenerate case. 
Such anomalies cannot be explained in the context of the degenerate = 2 model within 
the parameter space allowed by other neutrino oscillation experiments. However, we have 
argued that in the non-degenerate case, it is possible to obtain the pattern favoured by 



1 14], although the same tension between appearance and disappearance still remains. The 



parametrization of Casas- Ibarra [11] in the mini-seesaw region of the general model has 
been used to derive accurate approximations of the oscillation probabilities in this regime. 

The detailed constraints for the non degenerate case including CP violation will be 
considered in future work. We have not explored the restricted parameter space implied 
by a possible approximate lepton number symmetry, that can induce technically natural 
hierarchies in the spectrum. Such approximate symmetry could imply cancellations that 
could have been missed in our scan, and need to be searched for more carefully. Such 
possibility, as well as the model with increased level of complexity implied by an additional 
Weyl fermion, ur = 3, are interesting avenues to be explored. 
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